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Abstract. 

We present preliminary results from Smooth Particle Hydrodynamics (SPH) 
simulations of common envelope evolution. We qualitatively compare the in- 
teraction between a 0.9Mq red giant with two different companion masses: a 
0.05-Mq brown dwarf and a 0.25Mq white dwarf companion. 



1. Introduction 

Most stars (60%) are part of binary or multiple systems. If the binary separation 
is small enough, the more massive star will eventually engulf the companion as 
it expands during its giant phase, and form a so-called common envelope (CE) 
system. Friction then transfers orbital energy and angular momentum into the 
envelope, which may eventually lead to ejection of the entire envelope. Most 
tight binary systems are believed to have been in a CE phase at one point in 
their life: high- and low-mass X-ray binaries, double degenerate white dwarf and 
neutron star binaries, cataclysmic variables and supersoft X-ray sources. Many 
stellar types - such as hydrogen deficient stars, early R-stars, or Wolf-Rayet 
type central stars of planetary nebulae for example - have also been proposed 
to be the result of a CE phase, as is olated stellar evolution models are not abl e 
to fully explain their properties (e.g. lDe Marco et al.| [2003a; I zzard et al.l 12007). 



New studies have even propose d that the majority of p lanetary nebulae may be 



the result of a CE interaction (Moe De Marco 200 



The concep t of CE evolution has been around since the pioneering work of 
Paczvnskl ( 19761 ). Earlier simulation efforts ( Rasio &: Livid[l996l ; ISandquist et al. 



1998) were restricted to high mass companions due to the low numerical reso- 



lution achievable at that time. Only a handful low-resolution simulations are 
available as of now, covering only a tiny frac t ion o f the vast parameter space 
available. Only very recently, iRicker Sz Taam ( 20071 ) have restarted simulating 



these dynamic events at high resolution with the adaptive mesh refinement code 
FLASH, so far concentrating again on interactions with relatively high- mass 
companions. Theoretical considerations expect lower-mass companions (such 
as brown dwarfs, or even Jupiter-mass planets) to penetrate much deeper into 
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Figure 1 . Time sequence of SPH particle plots for a common envelope be- 
tween a 0.9Mq red giant and a 0.25-Mq white dwarf. The small white dots 
show individual SPH prarticles, while the grey filled circles show the degerate 
core of the red giant (center) and the companion (right, on the surface). The 
individual plots are spaced 40 days apart from to 200 days. Note how the 
companion transfers energy and angular momentum into the envelope, which 
is transported outward in a spiral density wave. At the end of this simulation, 
the binary has almost stalled its evolution, as there is only little material left 
at its radius to interact with, and most of the envelope corotates with the 
companion. 



the giant's envelope or even mer ge into the core, requirin g a muc h higher nu- 



merical resolution at the center, iDe Marco et al.1 (see e.g. l2003bl ). Yet, these 



binaries should be much more common and may easily prove powerful enough 
to shed the rather losely bound envelope of a red giant (RG) or asymptotic gi- 
ant branch (AGB) star. Here we will present a qualitative comparison between 
two CE evolution simulations involving a 0.9Mq red giant interacting with a 
low (0.05Mq brown dwarf, BD) and medium mass (0.25Mq white dwarf, WD) 
companion. The BD-RG system was chosen to be the progenitor system of a re- 
cently discovered tig ht binary system between a 0.05Mq BD and a 0.39Mq WD 
(|Maxted et al.ll2006l ). Our RG was chosen such that its degenerate core matches 
this WD mass. We will use this system as a validation case for our code, and 
thus expect the BD to spiral in as far as the observed 0.65i?Q separation. 



2. SPH Simulations 

We use the Smooth Particle Hydrodynamics (SPH) technique to simulate the CE 
evolution of a 0.9Mq RG interacting with two different companions: a medium 
mass companion (0.25Mq WD, Figure [1]) and a low mass companion (0.05Mq, 
Figure [2]). Both simulations start with the companion in a circular orbit directly 
on the surface of the non-corotating star at 83-Rq. We use 100k SPH particles 



Figure 2. Same as Figure[T] but for a 0.05Mq brown dwarf companion. The 
images are also spaced 40 days apart, but extend until 320 days of evolution. 
Note how the interaction is much gentler, as the companion can transfer less 
energy and angular momentum into the envelope. The companion sinks in 
much deeper, and is in fact still sinking rapidly at the end of the evolution. 
The most energetic phase of this interaction is yet to come. 



in these low-resolution te st runs, still a facto r of 2 higher than the last SPH 
simulations on this topic ( Rasio Livid fl996l ). 

Note how the higher mass WD companion is able to shed the outer layers 
of the star very quickly. Within only about 100 days, the binary orbit shrinks to 
a little more than a third (30-Rq) of the initial radius (Figure [31 left). However, 
after that, the evolution slows down significantly. Since part of the envelope has 
already been shed, there is less envelope material at this separation to transfer 
energy to. Additionally, the envelope slowly starts to corotate with the WD, 
making it even more difficult for the companion to interact. The lower-mass BD 
companion on the other hand interacts more gently with the envelope. Due to 
the lower orbital energy and angular momentum of the system, the companion 
cannot shed the outer layers as quickly. The BD slowly but steadily sinks deeper 
into the RG's envelope and is still sinking quickly at the current end of the 
simulation at around 350 days. We expect the BD to sink as deep as O.65.R0, as 
sugge sted by the recent observation of a similar BD-WD binary bv lMaxted et al 



(|2006r ). which we are trying to mimic here. Another major difference between 
the two simulations is the location of the ejected envelope material. The WD 
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Figure 3. Le/t: Separation between RG core and companion as a function 
of time for the for the 0.05Mq (grey) and 0.25Mq companion (black). Center 
and Right: Mass distribution as a function of time for the 0.05Mq (left) and 
0.25Mq companion (right). The area under the curve is proportional to the 
mass at that radius. Colors indicate the evolved time, ranging from (darkest) 
to 300days (lightest). The higher-mass companion expels the envelope much 
faster with the bulk of the mass already being at around 1000i?Q. For the 
lower-mass companion on the other hand, most of the mass still hangs around 
100.Rq. Material that is being ejected later at a higher speed has to "plow" 
through this wall first, potentially changing the dynamics of the evolution. 



companion interacts violently with the outer layers of the RG and quickly moves 
the bulk of the material beyond IOOORq (Figure EJ right). The BD on the 
other hand ejects the material more gently, and the bulk still "sits" just beyond 
100 Rq, just beyond the original surface of the RG (Figure El middle), and there 
is some evidence that some of it may fall back down on the star again. Thus we 
expected the subsequent evolution of both cases to behave quite differently, as 
the most energetic phase of the BD interaction is still to come as it sinks deeper 
into the potential well. Any material that will be ejected from then on will have 
to first "plow" through the wall of cooled material that has accumulated around 
the star. 

In a forthcoming paper, we will explore the dynamics of CE evolution with 
higher resolution SPH simulations, putting particular emphasis on the neglected 
parameter space for low-mass companions. 
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